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ABSTRACT 

We report the discovery of the first pulsating extremely low mass (ELM) white dwarf (WD), SDSS 
J184037.78+642312.3 (hereafter J1840). This DA (hydrogen-atmosphere) WD is by far the coolest 
and the lowest-mass pulsating WD, with T e ff=9100 ± 170 K and log g=6.22 ± 0.06, which corresponds 
to a mass ~ 0.17 Mq. This low-mass pulsating WD greatly extends the DAV (or ZZ Ceti) instability 
strip, effectively bridging the log g gap between WDs and main sequence stars. We detect high- 
amplitude variability in J1840 on timescales exceeding 4000 s, with a non-sinusoidal pulse shape. Our 
observations also suggest that the variability is multi-periodic. The star is in a 4.6 hr binary with 
another compact object, most likely another WD. Future, more extensive time-series photometry of 
this ELM WD offers the first opportunity to probe the interior of a low-mass, presumably He-core 
WD using the tools of asteroseismology. 

Subject headings: binaries: close - - Galaxy: stellar content - - Stars: individual: SDSS 
J184037.78+642312.3 — Stars: white dwarfs — variables: general 



1. INTRODUCTION 

Asteroseismology allows us to probe below the photo- 
sphere and into the interiors of stars. There are many 
pulsational instability strips on the Hertzsprung-Russell 
diagram, including the DAV (or ZZ Ceti) instability 
strip, driven by a hydrogen partial ionization zone in 
the hydrogen atmosphere (DA) WDs. Seismology us- 
ing the non-radial g-mode pulsations of DAVs enables 
us to constrain the mass, core and envelope composi- 
tion, rotation rate, and the behavior of convectio n in 
these objects (see reviews by I Winget fe~K epler 20Ql and 
iFontaine fc BrassardfeOOSl ) . 

The mass distribution of DA WDs in the SDSS shows a 
strong peak at 0.6 M^with ta ils toward higher and lower 
masses (|Tremblav et alJ[20Tlh : masses of individual WDs 
range from about 0.2 M to 1.3 M Q . The roughly 150 
DAVs known to date have masses 0.5 — 1.1 M Q , implying 
they all likely contain C/O-cores. Lower mass WDs are 
likely to pulsate as well. Howeve r, previous searches ha ve 
failed to detect such pulsations (|Steinfadt et al.ll2012l ). 

The galaxy is not old enough to produce low mass 
(< 0.5 M ) WDs through single-star evolution; these 
WDs are believed to be the product of binary evolu- 
tion. Indeed, radial velocity surveys of lo w-mass WDs 
indicate that most form in binary systems flMarsh et al.1 
119951 [Brown et aLll2011a» . Many short-period binaries 
go through one or two common-envelope phases, which 
may effectively remove enough mass to prevent ignition 
of He to C/O. However, there is little direct evidence 
that low-mass WDs have He-cores. But if they pulsate, 
as do their C/O-core brethren, we may differentiate their 
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interior structure. 

We have been engaged in an ongoing search for low- 
mass DAVs for many years at McDonald Observa- 
tory. The benefits of a search for a low-mass (and 
putatively He-core) D AV were recently emphasized by 
ISteinfadt et al.l (|2010l ). Should they pulsate in g- modes 
like the C/O-core DAVs, the eigenfunctions of ELM WDs 
would globally sample the interior, making the pulsations 
sensitive to core composition. Seismology may also allow 
us to constrain the hydrogen layer mass; this is vitally 
important since hydrogen burning is expected to be a 
major or even d ominant componen t of the luminosity of 
e.g..lPanei et al.ll2007t) . 



these stars 
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(|2012h 



outlined the null results of a 
search for pulsations in 12 low-mass WDs. We have 
extended a similar, systematic search for variable He- 
core WDs, armed with the many dozens of new ex- 
tremely low-mass (ELM, ~ 0.2 M&) WDs catalogued by 
the ELM Survey dBragn et al. l [2010 : iKilic et al.ll2011al: 
IBrown et ail 120121 iKilic et al.1 \20m . That search has 



yielded its first success. 

In this Letter, we report the discovery of pulsations 
in the ELM WD SDSS J184037.78+642312.3, and show 
that the photometric variations are likely multi-periodic. 
We have also included our null results for another eight 
low-mass WDs that were observed not to vary, to various 
detection limits. 

2. OBSERVATIONS 

IBrown et al.1 (|2012l ) present the spectroscopic discov- 
ery data for J1840 from the Blue Channel spectrograph 
on the 6.5m MMT. They use 37 separate spectra over 
more than a year to determine the system parameters. 
IBrown et a l. (2012) find that J1840 is in a 4.5912±0.0012 
hr (16528.32 ± 4.32 s) orbital period binary with a K — 
272 ±2 km s _1 radial velocity semi-amplitude. However, 
a significant alias exists at 3.85 hr. 

Model fits to the co-added spectra for J1840 yield T e g 
= 9140±170 K and log g = 6.16±0.06, which correspond 
to a mass of roughly 0.17 M@. Given the mass function 
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TABLE 1 

Frequency solutions for SDSS J1840+6423 



Fig. 1. — High-speed photometry of J1840 (in black) over three consecutive nights in 2011 October. The gap in the third night was 
caused by a passing cloud. The bottom red light curve, offset by —15%, shows the brightest comparison star in the field over the same 
period. Average 2a errors are shown offset with error bars. 

of the system, / = 0.399 ±0.009 M©, the minimum mass 
of the unseen companion is 0.64 M ; if the orbital in- 
clination is random, there is a 70% probability that the 
companion is a WD with < 1.4 Mq. We note, however, 
that the nature of the unseen companion has no bearing 
on the impact of this Letter. 

We obtained high speed photometric observations of 
J1840 at the McDonald Observatory over three consec- 
utive nights in 2011 October, for a total of more than 
5.5 hr of coverage. We used the Argos instrument, a 
frame-transfer CCD mounte d at the prime focus of th e 
2.1m Otto Struve telescope (jNather fc Mukadamir2004h . 
to obtain many 15 s exposures on this go = 18.8 mag 
WD. The seeing averaged 1.5" and transparency vari- 
ations were low, although our second and third nights 
were cut short by clouds. Observations were obtained 
through a 1mm BG40 filter to reduce sky noise. 

We performed weighted aperture photometry on the 
calibrated frames using the external IRAF package 
ccd-hsp written by Antonio Kanaan (t he reduction 
method is outlined in lKanaan et al.l 12002ft . We divided 
the sky-subtracted light curves by five brighter compar- 
ison stars in the field to allow for fluctuations in see- 
ing and cloud cover, and applied a timing correction to 
each observation to account for the motio n of the Earth 
around the barycenter of the solar system (|Stumpflll980t 
iThomoson fc Mullallvl[2009h . 

Figure [1] shows all 1365 Argos light curve points ob- 
tained for J1840 from 25 Oct 2011 to 27 Oct 2011. We 
also include the light curve of the brightest comparison 
star in the field, SDSS J184043. 21+642351.8, for refer- 
ence. 



Period 


Frequency 


Amplitude 


S/N 


00 


(MHz) 


(%) 




Multi-mode solution I 


4445.9 ± 1.4 


224.926 ± 0.070 


6.59 ± 0.19 


8.0 


2376 ± 57 


420 ± 10 


4.883 ± 0.83 


6.3 


1578.56 ± 0.37 


633.49 ± 0.15 


2.95 ± 0.25 


4.4 


Multi-mode solution II 


4445.3 ± 2.4 


224.96 ± 0.12 


7.6 ± 1.6 


9.4 


2376.07 ± 0.74 


420.86 ± 0.13 


4.817 ± 0.46 


6.3 


1578.70 ± 0.65 


633.43 ± 0.26 


2.831 ± 0.41 


4.3 


3930 ± 300 


254 ± 19 


2.7 ± 2.0 


3.4 


1164.15 ± 0.38 


859.00 ± 0.29 


1.78 ± 0.29 


3.3 


Single-mode harmonic solution 


4443.77 ± 0.80 


225.034 ± 0.041 


7.35 ± 0.19 


5.7 


2221.89 ± 0.40 


450.068 ± 0.081 


4.03 ± 0.19 


3.5 


1481.26 ± 0.27 


675.10 ± 0.12 


1.53 ± 0.18 


1.6 



3. ANALYSIS 

Our photometric data set is relatively short, as we 
caught J1840 just before it went behind the Sun. This 
limits the significance of the detected periods, and we 
eagerly anticipate further observations. Still, we have 
sufficient data to show convincingly that this low-mass 
WD is a multi-periodic variable star. 



The high-amplitude variability is easy to distinguish 
in the raw light curve (Figure [IJ , with more than 25% 
peak-to-peak variability. The highest peak in a Fourier 
transform (FT) of the brightest companion star in the 
field yields only a small signal (0.24% amplitude) at 
7340 ± 15 s, consistent with low-frequency noise from at- 
mospheric variability. Our apertures and sky annuli have 
been chosen to ensure there is no significant contamina- 
tion from the nearby star SDSS J184038. 73+642315.6, 
which is 7.0" away from our target. Thus the signal we 
are observing is intrinsic to the WD. Without evidence 
for a companion star or accretion from the spectra, the 
light curves, or existing broadband photometry, we are 
confident that the photometric variability results from 
pulsations on the surface of the ELM WD. 

The pulse shape appears non-sinusoidal, with a steep 
rise and decline. We first test whether a single mode 
(/i) and its harmonics (2/i + 3/i + ...) could reproduce 
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Fig. 2. — Fourier transforms of the light curve. The top panel 
shows the original data, the second panel the data after prewhiten- 
ing by the peak at 224.926 /iHz (4445.9 s), the third panel after 
prcwhitening by the peak at 420 /iHz (2376 s), and the fourth panel 
afterprewhitening by the peak at 633.49 /tHz (1578.56 s). See Ta- 
ble [TJ for our frequency solution. The dashed blue line shows 4 
times the noise level, as described in the text. The bottom panel 
shows the spectral window folded around 633.49 /iHz. 

the observed light curve. A nonlinear least squares fit 
with the highest peak in the FT and its fixed harmonics 
converges on 225.03 /iHz (4444 s) as the best parent mode 
(see the bottom of Tabled] for a full solution). 

However, our data from multiple nights rule out this 
scenario. The single-mode harmonic solution fits the two 
peaks during the first night well, but it fails to match the 
peak in our second night of data, predicting a maximum 
in the light curve more than 500 s too soon. 

A multi-mode solution (see the top panels of Table [lj 
where the frequencies are determined from the highest 
amplitudes in an FT of the entire dataset (see Figure [Uj 
improves the residuals by more than 20%. For more re- 
alistic estimates, the cited errors are not formal least- 
squares errors to the data but rather the product of 10 5 
Monte Carlo simulations o f perturbed data usi ng the 
software package Period04 (jLenz fc Bregerl l2005'l . The 
signal-to-noise calculation is based on the average ampli- 
tude of a 1000 /iHz box after pre-whitening by the three 
significant, highest-amplitude frequencies. 

A small but possibly significant amount of power re- 
mains around 254 and 859 /iHz after prcwhitening by 
the three significant periodicities (with S/N > 4), so we 
have included a frequency solution with these additional 
periods in Table [TJ This set of period solutions is by no 
means exhaustive, but it establishes that this ELM WD 
is variable and multi-periodic. 



4. DISCUSSION 

We announce the discovery of the first pulsating ex- 
tremely low mass white dwarf, J 1840, which is both the 
coolest and lowest-mass WD known to pulsate. The ob- 
ject offers, for the first time, an opportunity to explore 
the interior of a putative He-core WD using asteroseis- 
mology. 

Asteroseismology of ELM WDs will help constrain the 
thickness of the surface hydrogen layers in these low- 
mass WDs. There are several millisecond pulsars with 
ELM WD companions, and the cooling ages of such 
companions can be used to calibrate the spin-down ages 
of these pulsars. However, current evolutionary mod- 
els for ELM WDs are relatively unconstrained. For 
masses M > 0.17 A/©, diffusion-induced hydrogen-shell 
flashes take place, which yield small hydrogen envelope s 
(|Althaus et al.|[200lt iPanei et al.ll2007t iKilic et al.lf2010h . 
The models with M < 0.17 M do not experience ther- 
monuclear flashes. As a result, they have massive hydro- 
gen envelopes, larger radii, lower surface gravities, and 
they are predicted to evolve much more slowly compared 
to more massive WDs. If enough modes are excited to 
observability, we hope to directly constrain the hydrogen 
layer mass. 

A plethora of excited modes would also allow for mea- 
suring the mean period spacing, which is a sensitive func- 
tion of the mass of the star, and is also sli ghtly dependent 
on the core mass fraction. The models of Stcinfad t et al.l 
(pOlTl found a mean period spacing of ~ 89 s for £ = 1 
(/-modes of a 0.17 M ELM WD, about a factor of two 
larger than the observed 47 ± 12 s period spacing for 
£ = 1 q-modes in the c ool C/O-core counterpart G29-38 
(jKleinman et al.lll998| ). 

We do expect this non-radial pulsator to be multi- 
periodic even if just one pulsation mode and its non- 
linear harmonics are amplified to observability: J 1840 
is in a relatively close binary with an unseen compact 
object. This companion will influence the light curve 
in many ways. Although we lack the sensitivity and 
full phase coverage needed for a detection, we expect 
a ~ 0.3% amplitude Doppler beaming signal at the or- 
bital period given the effective temperature of J184 and 
the radial- velocity amplitude ([Shporer et alj |2010). The 
companion will also induce tidal distortions on the pri- 
mary, as seen m many ot her ELM WDs (IKilic et al.1 
120 lit* iBrown et al.l IHermes et all 12012ft . How- 

ever, assunfiirgtha^ J L1840 has a radius of roughly 0.054 
Rq (|Panei et al.l 120071) . ellipsoidal variations should be 
small er than 0.1% for even th e highest possible inclina- 
tions (jMorris fc Naftilanl[T99l . 

The companion's effect on the rotation period is per- 
haps more s ignificant. The ~ 1.7 Gyr cooling age of this 
ELM WD (|Panei et all 120071) may be longer than the 
synch ronization timesca le for such a short-period binary 
(e.g., iClaret et al1ll995D . If so, and the ELM WD is ro- 
tating at the orbital period of 4.6 hr, then non-radial pul- 
sations will be subsequently subject to rotational split- 
tings determined by their modal degree. For example, 
if the 225 /iHz mode is an £ = 1 mode, a 4.6 hr rota- 
tion rate would cause it to b e split by about 30 /iHz, 
assuming solid-body rotation (jUnno et al.lll989T) . How- 
ever, synchronization is not a guarantee: Recent Kepler 
observations of a close sdB+dM binary found the pri- 
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Fig. 3. — Search parameter space for pulsations in low-mass 
WDs. Known DAVs are included as purple dots, and our new 
pulsator is labeled and marked as a magenta dot. DA WDs ob- 
serve d not to vary to at l east 1% are marked with an X; those 
from Stcinfadt ct al. (2012) are marked in gray and those from our 
search (see Table [2{ in black. To guide the eye, we have marked 
the rough boundaries of the observed DAV instability strips as 
dotted line s; the short, blu e dotte d line marks the blue-edge of the 
theoretical Stcinfadt ct al. (2010) instability strip. 

mary r otating much slow er than the ~ 9.5 hr orbital 
period (|Pablo et al.ll2012l ). An analysis of the rotational 
splittings in the pulsations of J1840 will test tidal syn- 
chronization in this system, thereby probing the rigidness 
of this ELM WD. 

We also hope to use the nonlinearities in this non- 
sinusoidal light curve to constrai n the size of the 
convection zone of this WD (e.g., Montgo mery! 120051 : 
iMontgomerv et al.l l2010h . Since ELM WDs exist in a 
completely new regime of parameter space than do typ- 
ical C/O core WDs, this analysis will provide important 
independent constraints on the structure of their outer 
layers, as well as providing a measure of the convective 
efficiency in this new regime. 

We discovered the pulsations in J 1840 as part of a sys- 
tematic search for pulsatio ns in low-mass WDs, recentl y 
energized by the work of iSteinfadt et al.l (|2010l I2012I ). 
We have explored a vast area of the log g — T c g parame- 
ter space for variability (see Figure[3]), and compiled a list 
of our first eight null results in Table [2] These null obser- 
vations were reduced and analyzed in an identical man- 
ner to those outlined in Section [2j and the atmospheric 
parameters and their formal errors were determined by 
the references cited. 

Our newfound pulsator J 1840 occupies a new, cooler 
region of the DAV instability strip, which may be an ex- 
tension of the C/O-core region. The highest-amplitude 



period observed (> 4400 s) is the longest, to date, of any 
period observed in a DAV (|Mukadam et alj 12006). This 
makes sense qualitatively, as we would expect the periods 
of pulsation modes to roughly scale with the dynamical 
timesc ale for the whole star, P oc p -1 / 2 . 0st ensen et al.l 
(2010) find longer pulsation periods in lower tempera- 
ture and surface gravity subdwarf B stars. A similar 
trend would explain the relatively long pulsation period 
of J 1840. The pe riod observed is a fac tor of 4 longer 
than the models of Stcinfad t et al.l (|2010l ) predict for an 
£ = 1, k = 10 mode of a 0.17 M ELM WD; in the con- 
text of these models, this mode likely has an extremely 
high radial overtone, higher than that of any dominant 
mode observed in a normal-mass DAV. 

We note, however, that the driving mechanism for pul- 
sation is largely unknown for this star. It is natural to 
assume that the same mechanism of convective driving 
that operates in the C/O-cor e DAVs (|Brickhilllll99lH Wi] 
119981 : iGoldreich k, Wul 11999! ) is also responsible for the 
pulsations of low-mass WDs. This mechanism is based 
on the assumption that the convective turnover timescale 
for a fluid element, t to , is much smaller than the oscilla- 
tion periods, Pj. For the logg ~ 8 DAVs, t to ~ 0.1-1 sec 
and Pi ~ 100-1000 sec, so t to -C Pi is satisfied. 

To estimate how this timescale scales with g we note 
that tt ~ l/v, where / is the mixing length and v is the 
velocity of the convective fluid elements. Employing a 
simplified version of mixing length theory (ignoring ra- 
diative losses) we find that v oc (gF/ p) 1 ^ 3 , where F is the 
stellar flux and p is the mass d ensity at the b ase of the 
star's surface convection zone (Mihalas 1978). Taking 
the mixing length as a factor of order one times a pres- 
sure scale height, we find I ~ c 2 s /g oc T/g. Putting these 
results together we find t to oc T(p/F) 1 / 3 g- 4 / 3 - g- 4 / 3 . 
Thus, a logg ~ 6 object would have t to ~ 50-500 sec. 
In order for convective driving to operate, p ^> t to is 
required. Perhaps this is at least a partial explanation 
for the very long period (> 4000 s) seen in this pulsator. 
It may also set a lower limit to periods in this DAV of 
~ 500-1000 s, although more detailed models will be 
needed to confirm this. Another potential source of driv- 
ing is the e mechanism, i.e., driving due to the modula- 
tion of H burning in the envelope. We will address this 
additional mechanism after more extensive observations. 

We look forward to a coordinated effort for extensive 
follow-up observations and analysis of this exciting new 
pulsating WD, and to the discovery of additional pulsat- 
ing ELM WDs in order to better understand this new 
(or extended) DAV instability strip. 
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